home *** CD-ROM | disk | FTP | other *** search
/ AmigActive 24 / AACD 24.iso / AACD / Sound / LAME / WarpOS / src / mpglib / tabinit.c < prev    next >
Encoding:
C/C++ Source or Header  |  2001-02-05  |  4.6 KB  |  124 lines
  1. /* $Id: tabinit.c,v 1.11 2001/01/05 15:20:34 aleidinger Exp $ */
  2.  
  3. #ifdef HAVE_CONFIG_H
  4. # include <config.h>
  5. #endif
  6.  
  7. #include <stdlib.h>
  8. #include "tabinit.h"
  9. #include "mpg123.h"
  10.  
  11. #ifdef WITH_DMALLOC
  12. #include <dmalloc.h>
  13. #endif
  14.  
  15. real decwin[512+32];
  16. static real cos64[16],cos32[8],cos16[4],cos8[2],cos4[1];
  17. real *pnts[] = { cos64,cos32,cos16,cos8,cos4 };
  18.  
  19. const double dewin[512] = {
  20.    0.000000000,-0.000015259,-0.000015259,-0.000015259,
  21.   -0.000015259,-0.000015259,-0.000015259,-0.000030518,
  22.   -0.000030518,-0.000030518,-0.000030518,-0.000045776,
  23.   -0.000045776,-0.000061035,-0.000061035,-0.000076294,
  24.   -0.000076294,-0.000091553,-0.000106812,-0.000106812,
  25.   -0.000122070,-0.000137329,-0.000152588,-0.000167847,
  26.   -0.000198364,-0.000213623,-0.000244141,-0.000259399,
  27.   -0.000289917,-0.000320435,-0.000366211,-0.000396729,
  28.   -0.000442505,-0.000473022,-0.000534058,-0.000579834,
  29.   -0.000625610,-0.000686646,-0.000747681,-0.000808716,
  30.   -0.000885010,-0.000961304,-0.001037598,-0.001113892,
  31.   -0.001205444,-0.001296997,-0.001388550,-0.001480103,
  32.   -0.001586914,-0.001693726,-0.001785278,-0.001907349,
  33.   -0.002014160,-0.002120972,-0.002243042,-0.002349854,
  34.   -0.002456665,-0.002578735,-0.002685547,-0.002792358,
  35.   -0.002899170,-0.002990723,-0.003082275,-0.003173828,
  36.   -0.003250122,-0.003326416,-0.003387451,-0.003433228,
  37.   -0.003463745,-0.003479004,-0.003479004,-0.003463745,
  38.   -0.003417969,-0.003372192,-0.003280640,-0.003173828,
  39.   -0.003051758,-0.002883911,-0.002700806,-0.002487183,
  40.   -0.002227783,-0.001937866,-0.001617432,-0.001266479,
  41.   -0.000869751,-0.000442505, 0.000030518, 0.000549316,
  42.    0.001098633, 0.001693726, 0.002334595, 0.003005981,
  43.    0.003723145, 0.004486084, 0.005294800, 0.006118774,
  44.    0.007003784, 0.007919312, 0.008865356, 0.009841919,
  45.    0.010848999, 0.011886597, 0.012939453, 0.014022827,
  46.    0.015121460, 0.016235352, 0.017349243, 0.018463135,
  47.    0.019577026, 0.020690918, 0.021789551, 0.022857666,
  48.    0.023910522, 0.024932861, 0.025909424, 0.026840210,
  49.    0.027725220, 0.028533936, 0.029281616, 0.029937744,
  50.    0.030532837, 0.031005859, 0.031387329, 0.031661987,
  51.    0.031814575, 0.031845093, 0.031738281, 0.031478882,
  52.    0.031082153, 0.030517578, 0.029785156, 0.028884888,
  53.    0.027801514, 0.026535034, 0.025085449, 0.023422241,
  54.    0.021575928, 0.019531250, 0.017257690, 0.014801025,
  55.    0.012115479, 0.009231567, 0.006134033, 0.002822876,
  56.   -0.000686646,-0.004394531,-0.008316040,-0.012420654,
  57.   -0.016708374,-0.021179199,-0.025817871,-0.030609131,
  58.   -0.035552979,-0.040634155,-0.045837402,-0.051132202,
  59.   -0.056533813,-0.061996460,-0.067520142,-0.073059082,
  60.   -0.078628540,-0.084182739,-0.089706421,-0.095169067,
  61.   -0.100540161,-0.105819702,-0.110946655,-0.115921021,
  62.   -0.120697021,-0.125259399,-0.129562378,-0.133590698,
  63.   -0.137298584,-0.140670776,-0.143676758,-0.146255493,
  64.   -0.148422241,-0.150115967,-0.151306152,-0.151962280,
  65.   -0.152069092,-0.151596069,-0.150497437,-0.148773193,
  66.   -0.146362305,-0.143264771,-0.139450073,-0.134887695,
  67.   -0.129577637,-0.123474121,-0.116577148,-0.108856201,
  68.   -0.100311279,-0.090927124,-0.080688477,-0.069595337,
  69.   -0.057617187,-0.044784546,-0.031082153,-0.016510010,
  70.   -0.001068115, 0.015228271, 0.032379150, 0.050354004,
  71.    0.069168091, 0.088775635, 0.109161377, 0.130310059,
  72.    0.152206421, 0.174789429, 0.198059082, 0.221984863,
  73.    0.246505737, 0.271591187, 0.297210693, 0.323318481,
  74.    0.349868774, 0.376800537, 0.404083252, 0.431655884,
  75.    0.459472656, 0.487472534, 0.515609741, 0.543823242,
  76.    0.572036743, 0.600219727, 0.628295898, 0.656219482,
  77.    0.683914185, 0.711318970, 0.738372803, 0.765029907,
  78.    0.791213989, 0.816864014, 0.841949463, 0.866363525,
  79.    0.890090942, 0.913055420, 0.935195923, 0.956481934,
  80.    0.976852417, 0.996246338, 1.014617920, 1.031936646,
  81.    1.048156738, 1.063217163, 1.077117920, 1.089782715,
  82.    1.101211548, 1.111373901, 1.120223999, 1.127746582,
  83.    1.133926392, 1.138763428, 1.142211914, 1.144287109,
  84.    1.144989014
  85. };
  86.  
  87. void make_decode_tables(long scaleval)
  88. {
  89.   int i,j,k,kr,divv;
  90.   real *table,*costab;
  91.  
  92.   
  93.   for(i=0;i<5;i++)
  94.   {
  95.     kr=0x10>>i; divv=0x40>>i;
  96.     costab = pnts[i];
  97.     for(k=0;k<kr;k++)
  98.       costab[k] = (real)( 1.0 / (2.0 * cos(M_PI * ((double) k * 2.0 + 1.0) / (double) divv)) );
  99.   }
  100.  
  101.   table = decwin;
  102.   scaleval = -scaleval;
  103.   for(i=0,j=0;i<256;i++,j++,table+=32)
  104.   {
  105.     if(table < decwin+512+16)
  106.       table[16] = table[0] = (real)( dewin[j] * scaleval );
  107.     if(i % 32 == 31)
  108.       table -= 1023;
  109.     if(i % 64 == 63)
  110.       scaleval = - scaleval;
  111.   }
  112.  
  113.   for( /* i=256 */ ;i<512;i++,j--,table+=32)
  114.   {
  115.     if(table < decwin+512+16)
  116.       table[16] = table[0] = (real) ( dewin[j] * scaleval );
  117.     if(i % 32 == 31)
  118.       table -= 1023;
  119.     if(i % 64 == 63)
  120.       scaleval = - scaleval;
  121.   }
  122. }
  123.  
  124.